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DIFFERENTIAL EQUATI0N8 IK AIRPLANE MECHANICS.* 

By 

li. T. Carleman. 

For determining the motion of an airplane, we have adopted 
the hypothesis that the reactions of the air depend entirely on 
the relative speed of the airplane. Even if we adopt the sim- 
plest laws of resistance, we obtain differential equations whioh 
we can not integrate explicitly. If we oonfine ourselves to the 
motion in a vertical plane and, at ths same tine, assume a con- 
stant angle of attaok, we still obtain differential equations 
which oan not be integrated by elementary methods. 

In the following paragraphs, we Till first draw sor.e con- 
clusions of purely theoretioal interest, from the general equa- 
tions of motion. At the end, we will oonsider the motion of an 
airplane, with the engine dead and with the assumption that the 
angle of attach remains oonstant. Thus we arrive at a simple 
result, whioh osn be rendered praotically utilizable for deter- 
mining the trajectory of an airplane desoending at a oonstant 
stesring angle. 

Let us assume that the airplane moves in its plane of sycrs- 
try, considered vertical. Let x and y represent a system of 
coordinates in this plane, the axis of the x lines being hori- 
zontal and the axle of the y lines being vertloal and upward. 

Let vs designate by u and t the projeotions of the velocity w 
* From "La Teohnique Aeronautique, 11 ' Hay, 1931. 
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of the oenter of gravity of the airplane on the axes of x and y 
and by the angle \jr the angle of the velooity with the axis of 
the x lines. Leu us designate by If the mass of the airplane 
end by IE its moment of inertia with relation to an axis pass- 
ing through the oenter of gravity at right angles to the plane 
xy. The foroes acting on the airplane are'* 
Its weight Hg direoted toward -y: 

Propeller thrust T . of whioh we will designate by Ha and 
KB the components in the direction of w and at right angles 
to w; 

Drag, Ear^ in the direction -wj 
Lift. M b w a , at right angles to w. 

Under these conditions the theory of the notion of the oen- 
ter of gravity giveb vis 

M§g = - K a w 3 oos ■}■ M A oos ^ - II b w a Bin if - Iff B sin if. 

m 

= - M a w a sin >|r + K A sin if + Iff b w 3 oos i|r +11 B cos if • 

-Mg dt . 

By substituting 

cos if = |J sin if = J 

re have 

. (1) = - w (au+bv) 

(3) ^ = - w (a v - "b u) + Ay + ^ - g 
w = y-ua + v a . 
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If we project these forces on the tangent and on the per- 
pendicular to the trajectory, ^& obtain 

(3) $f - a-wa - g cin * + A. 

(4) = bw* - g oos i|r + B. 

By letting 

a = p sin a A = K oos p 
b = p 003 a B = K 3in 0 

T7e can, by introducing imaginary quantifies, combine equations 
(l) and (2) into a single equation 

( 5 ) flfaj It ) = (i w p 8 ia + *£?)(u + iv) - gi. 

The moment of the air resistances about the oenter of grav- 
ity can be expressed in the form 

1! C-( 9) - II H(«p). 

in nhioh G and H are the psriodio funotions of the angle of at- 

•taok <p, chioh varies with the steering angle. l?e then have, 
according to the theory of motion about the oenter of gravity, 

(6) I * 8 CjLfrl ) .= 0(9) - w H(q>). 

' Let us assume that the influence of altitude may be neglect- 
ed. By multiplying equation (l) by u and equation (2) by v 
and adding the results, we obtain 



- 4 - 

(?) -i = - aw 3 + Aw - gv. 

The quantity a, neoessarily implying a lower positive limi J 
a' , is obtained from equation (7). 



( 8 ) 1 flEf < - w / a 'w a - (P + g) "\ 
■ 3 dt V j 



HP representing the maxinum propeller thrust, from whioh we rea- 
son that the speed oan not go on increasing Indefinitely. To be 
more exact, it must remain below 



/ 



a' 



in whioh a" is any quantity smaller than a*. 

By means of equation (l) we can easily demonstrate the fol- 
lowing theorem, which is practically self-evident. If the engine 
is dead, it is not -possible fry R.jry ma neuver of the elevator to 
fre ?p th e ?irelane above anv fixed horizontal line for anv arbitra- 
ry period of £laa. T7e find, in- fact, by integrating equation (7; 
batween t 0 and t(A and B being zero) 

•i (w t s - w Q 3 ) = - /*■ aw 3 dt - g (y t - y Q ). 
a t Q 

Taking into aooount the inequality a > a' > 0, the hypothesis 
y-fc > 0 gives ua 

t T 
a' ^ w 3 dt < ^ w 0 3 + gy 0 . 

It follows that the integral 

f° w 3 dt 
t o 
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is convergent. Remarking that by reason of equations (7) and 

(p) dEi i s limited, we reason from this that' w tends towards 
dt 

aero, when t tends toward infinity. Henoe we have also 
(fe) lim 7 = 0, lim u = 0 

+ = V> t = » 

Consequently, there is an infinity of values of t, tj , % . . . t n 
such that we have 

(10) lim /dv\ m Q lim /du\ = 0j lim t a =<x> 
r. =«' dt / tn n = a Ut t = o» 

Py subeti-Gutlng t = t n in the equation 

4? = - w (av - bu) - g 

dt 

and by raking n tend toward infinity, we arrive at the absurdi- 
ty g - 0. The hypothesis y^ > 0 is therefore inadmissible. 

Let us now assume the thrust and angle of attaok to remain 
constant. On dividing equation (6) by r a , we see that for a 
constant steering angle of the elevator, this hypothesis con- 
cerning the angle of at-saok is judtlfied in proportion as I. 
and K ore smaller and w is larger. From equations (l) and 
(3) we deduoe the differential equation 

(id *>- 



(u 3 + v 3 ) (au + bv) - An + Bv 

dv_ 

(u s + v 8 ) (av - bu) - Av - Bu + g /u a + v a 



f 
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in which a, b, A and B axe oon3tantB. Once this equation is 
integrated, t, x and y are obtained "by quadratures. !Te can 
accordingly confine ourselves to the study of .the integral curves 
of equation (ll), that is, to determining the hodo graph of the 
motion. Ve oan apply to equation (ll) the methods given by l£r. 
Poinoare* in his memoirs : "On ourves defined by a differential 
equation" (Journal de Ifathematlques pures et appllquees, ie81- 
1886). tte see that equations (l) and (3) define u and v for 
every value of t, -when- the initial velocity UoV 0 is given. 
On oonsidering u and r as cartesian coordinates of a point 
in the plane, the point (u,v) desoribes a certain characteris- 
tic curve . -hen t varies from 0 to infinity. T?e pall "sin- 
gular- points" those Trhere the- second members of equations (l) 
and (s) disappear at the same time or become discontinuous. In 
the case under consideration, there is, eside from the origin, 
only one other singular point of finite distance. Tfe have al- 
ready seen (inequality 8) that the characteristics remain at a 
finite distanoe. Under this condition, the theory of Poinoare' 
shows us that a knowledge of the singular points and olosed char- 
acteristics (limited oyoles) at a finite distanoe suffices for 
finding the oourse of the integral ourves, when t tends toxrard 
infinity. 

Instead of going deeply into this study here, ve shall take 
up by a simpler method, the oase where the propeller thrust is 
zero. Equation (ll) is reduced to 
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/_ a x dtt _ dv 

xt (au + fcv) w (av - ba) + g 

For 0=0,* that is to say. a = 0, b « p, this equation is 
readily integrated by means of the elementary functions, tfe have, 
in faot, 

du „ dv xr d w 

p vw - p uw + g gv 

Hence 

gdu = p w s dw 
p ^ - gu = C (0- arbitrary oonstant) 

By varying 0 we obtain a system S of olosed, non-interseoting 
cmrves,. which embrace the whcle plane. 

In the case a = tt/3, corresponding to the motion of a sphere 
in a resisting medium, equation (13) oan be integrated by quadra- 
tures, as demonstrated by Liouville. 

Although we do not know for a any general integral of equa- 
tion (13), we oan construct its integral curves by a simple method. 
?& have, in faot, the following theorem. The integral curves of 
equation (13) are the trajectories of the angle a of the system of 
S ourves. 

p (u 3 -Mr 3 ) 3 °- gu = constant 

provided we consider these trajectories with referenoe to the reo- 

* In his article "Le vol Ae'rien 1 ' (Aerial Flight) , llr. Lanohester 
disouosed the fall of an airplane under this hypothesis, whioh re- 
turns to the assumption that the airplane does not lose energy. 



t angular axes which make the angle a (in the positive direction) 
with the original trajectories (see figure). 

For the demonstration, we write the differential equations in 
the following condensed form 

■ v (13) d (tt t ft ) = iw W fe (u + iv) - gi 

whioh is obtained from equation (5) by letting E *= 0. By making 
a = o in equation (13) we find, for the S system, the differen- 
tial equation 

(14) fl (u . + It) = iw p (u + iv) - gi. 

dt 

tfs obtain the differential equation of the trajectories of the an- 
gle a by multiplying the second member of equation (14) by e . 

Henoe 

( 15 ) fl fa d + 3a ) = iw p e ia (u + iv) - gie iQ 

The change of axes indioated in the above statement is obtained by 
multiplying u + iv in equation (15) by e ict . By dividing the 
ratio thus obtained by e ia , we return to equation (13) , whioh dem- 
onstrates the theorem. 

All the characteristics tend toward the point 



u = /■£ oos tt, v = - /§L sin -a 
J P J P 



whioh corresponds to a reotilinear and uniform descent.* If the ln- 
* If an integral curve passes near the point u = v = 0, we can 
hardly expect that it will correspond to the real motion of the air- 
plane, because the hypothesis that the angle of attaok remains con- 
stant is not justified for low speeds. 



tegral curve cuts the axis of the negative u's, n times, the air- 
plane executes n loops. We come, therefore, to the following 
conclusion. Whatever be the initial velocity, the airplane, af-csr 
executing, if necessary t a finite number of loops, acquires a mo- 
tion which approaches indefinitely a state of rectilinear and uni- 
form descent. 

Translated by National Advisory Committee for Aeronautios. 




